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Abstract. XMM-Newton observations of the outskirts of the Coma cluster of galaxies confirm the existence of a 
soft X-ray excess claimed previously and show it comes from warm thermal emission. Our data provide a robust 
estimate of its temperature (~ 0.2 keV) and oxygen abundance (~ 0.1 solar). Using a combination of XMM- 
Newton and ROSAT All-Sky Survey data, we rule out a Galactic origin of the soft X-ray emission. Associating 
this emission with a 20 Mpc region in front of Coma, seen in the skewness of its galaxy velocity distribution, yields 
an estimate of the density of the warm gas of ~ 50/baryonPcriticai, where /baryon is the baryon fraction of the gas 
and pcriticai is the critical density needed to halt the expansion of the universe. Our measurement of the gas mass 
associated with the warm emission strongly support its nonvirialized nature, suggesting that we are observing 
the warm- hot intergalactic medium (WHIM). Our measurements provide a direct estimate of the O, Ne and Fe 
abundance of the WHIM. Differences with the reported Ne/0 ratio for some OVI absorbers hints at a different 
origin of the OVI absorbers and the Coma filament. We argue that the Coma filament has likely been preheated, 
but at a substantially lower level compared to what is seen in the outskirts of groups. The thermodynamic state 
of the gas in the Coma filament reduces the star-formation rate in the embedded spiral galaxies, providing an 
explanation for the presence of passive spirals observed in this and other clusters. 

Key words, galaxies: clusters; individual: Coma — cosmology: observations — intergalactic medium — large-scale 
structure of universe 



1. Introduction 

The total amount of baryons directly observed in the local 
' Universe is only 35% of the value implied from primor- 
. dial nucleosynthesis, with that observed at high-redshift 
' (Fukugita et al. 1998, Table 3 excluding line 7') and that 
■ observed within local closed box systems, such as clusters 
' of galaxies (Briel et al. 1992). Numerical hydrodynamic 
simulations (e.g. Cen & Ostriker, 1999; Dave et al., 2001) 
indicate that about half of these missing baryons reside in 
a warm- hot intergalactic medium (WHIM). The tempera- 
ture of the WHIM is expected to be between 10^ and 10*" 
K with an overdensity with respect to the mean of (5 ^ 20. 
Most of the WHIM baryons are expected to reside outside 
virialized regions in diffuse large-scale filaments. 

Several surveys have been undertaken to find the 
WHIM in absorption. Tripp et al. (2000) have claimed 
a substantial fraction of the missing baryons in the form 
of OVI absorbers. However, their observations lack infor- 
mation on the ionization equilibrium as well as the O 
abundance required to link the abundance of OVI ions 
to underlying O mass and then to the amount of baryons. 
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Chandra and XMM grating results on detection of OVII 
and OVIII absorption lines (Nicastro et al. 2002; Fang et 
al. 2002; Mathur et al. 2003; Fang et al. 2003; Rasmussen 
et al. 2003) allowed for the first time an estimate of the 
ionization state of the gas, suggesting that the OVI ab- 
sorbers account for over 80% of the local baryons. Still, 
such claims rely on the assumed O abundance. 

Searches for the WHIM in emission with ROSAT have 
found some association of soft X-rays with galactic fila- 
ments (Wang et al. 1997; KuU & Bohringer 1999; Scharf et 
al. 2000; Zappacosta et al. 2002). Upper limits have been 
placed on emission by material connecting rich clusters of 
galaxies (Briel & Henry 1995) using stacked analysis of 
ROSAT All-Sky Survey data. The corresponding limit on 
the emitting gas density is S < 375, assuming a tempera- 
ture of 1 keV and an iron abundance of 0.3 solar. These 
observations may also be used to constrain the WHIM be- 
cause most of it is expected to be in filaments connecting 
clusters. 

The low filament densities implied by the weak absorp- 
tion and soft X-ray emission, and expected from typical 
overdensities in simulations, complicates WHIM detection 
and study, unless the filament containing it is viewed along 
the line of sight. A combination of both the requirement 
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Filament overdensity relative to the critical, (p/pcrit — 1), 
20 Mpc, area of the extraction regions in Fig. and rie = 8 
" Assuming a solar abundance ratio, 
flux, 10~^ photons cm~^ s~^ arcmin"'^ 
that a filament be so viewed and the expectation that 
clusters are at the intersection of filaments, indicates that 
searches for WHIM near clusters of galaxies may prove 
fruitful. The soft X-ray excess above the harder intraclus- 
ter medium (ICM) emission reported for some clusters 
(Lieu et al. 1996; Bonamente et al. 2002) may be signal- 
ing a fortunate orientation of a filament containing WHIM 
at those clusters. 




Fig. 1. An image of the Coma cluster in the 0.5-2 keV band, 
showing the positions of the spectral extraction regions (solid 
ellipses) with names according to the XMM-Newton survey 
notation. The coordinate grid marks R.A., Dec. (J2000.0). 

In this Paper we carry out an X-ray spectroscopic 
study of the outskirts of the Coma cluster, which has a 
soft excess, in order to search for and characterize any 
WHIM at this location. Our observations are from the mo- 
saic made by XMM-Newton. Kaastra et al. (2003) report 
a similar study of the center of the Coma cluster also us- 
ing data from the XMM mosaic. Not knowing the relative 
spatial distributions of the WHIM and ICM, it is difficult 
to predict whether observing the center or outskirts of a 
cluster will maximize the contrast of any warm compo- 
nent relative to the dominant hot component. However 
both studies do detect a warm component. We adopt a 
Coma cluster redshift of 0.023, iJo = 70 km s'^ Mpc-\ 



estimated by assuming a 90 Mpc distance, projected length of 
X lO"'^ cm^'^ in correspondence to pcrit for h=0.7 

and quote errorbars at the 68% confidence level. One de- 
gree corresponds to 1.67 Mpc. 

2. Observations 

The initial results of the XMM-Newton (Jansen et al. 
2001) performance verification observations of the Coma 
cluster are reported in Briel et al. (2001), Arnaud et al. 
(2001) and Neumann et al. (2001). In addition to the ob- 
servations reported there, three additional observations 
have been carried out, completing the planned survey of 
the Coma cluster. We have renamed the Coma-bkgd field 
from Briel et al. (2001) Coma-O in order to avoid confu- 
sion with the additional background pointings required to 
analyze these data. 

Since the Coma cluster is a large diffuse source, esti- 
mating the background is non-trivial. In general we fol- 
low the method of Lumb et al. (2002). The observations 
reported here are performed with the EPIC pn detector 
(Striider et al. 2001) using its medium filter. This config- 
uration leads to a different detector count rate from the 
cosmic (CXB) and Milky Way X-ray backgrounds in the 
soft band, compared to the usual thin filter. The Milky 
Way background is further composed of the Local Bubble 
and halo components. We used a 100 ksec observation of 
the point source APM08279+5255 (PI G. Hasinger), per- 
formed with the same (medium) filter, to measure these 
background components and a filter wheel closed observa- 
tion to measure the internal or particle background. We 
used the APEC plasma code (Smith et al. 2001) to fit 
the Milky Way background components and a power-law 
(r = 1.45) to fit the CXB. The element abundance in the 
APEC code was fixed to solar and temperature was set to 
0.07 keV for the Local Bubble and left as a free parameter 
for the halo, yielding a 0.22 keV best-fit value. To ensure 
the compatibility of the CXB, we excised from the spectral 
extraction region the APM quasar and the X-ray sources 
identified with Coma galaxies (which otherwise increase 
the apparent CXB flux by 10%). Since there is a mis- 
match in the Nh toward the Coma and the APM flelds, 
we use the same spectral shape and normalization of the 
halo and CXB components, but changed their from 
3.9±0.3xl02° cm-2 (APM) togxlO^^ cm^^ (Coma). The 
Local Bubble is not subject to this absorption. We also 
checked that the derived results on the soft X-ray emis- 
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sion from Coma do not depend on the possible large-scale 
variation in the intensity of the Local Bubble component. 

For the spectral analysis we use the 0.4-6 keV range 
and for the estimate of the detector background the 9.5- 
15 keV range. We do not use the energies below 0.4 keV, 
to avoid a complicated treatment of detector background 
associated with the electronic noise (Lumb et al. 2002). 
This precludes us from studying the very soft emission. 
Although, we could potentially use the energies up to 10 
keV in the analysis, we are able to sufficiently constrain 
the hot component of the Coma cluster already without 
using these energy bands, where removal of the detector 
background becomes critical. 

The vignetting correction is performed taking the 
source extent and a recent pn telescope (XRT3) vignetting 
calibration (Lumb et al. 2003) into account, which is 
mostly important for the absolute flux determination, 
given our choice of the energy range. Systematic uncer- 
tainty of the flux determination is 4% for both pn and 
MOS. 

In addition to a thorough evaluation of the back- 
ground, we further reduce the systematics of our study 
by using the results of recent calibration efforts of the pn 
team of MPE on soft energy response of the pn detector 
(Haberl et al. 2002, XMMSAS-5.4 release). Furthermore, 
to improve on the detection statistics of the warm emis- 
sion and to reduce the systematic effects of subtraction 
of the cluster emission, we concentrate on the outskirts of 
the Coma cluster. 

We investigated some systematic errors in our analysis 
by determining the effects of using experimental software 
or non-standard values of some parameters. We do not use 
these modifications in our analysis, however. The effect 
of pn energy scale calibration was evaluated by compar- 
ing the results of the standard analysis with experimental 
software that corrects the zero point of the energy scale for 
each pixel of the pn. We found some subtle differences, e.g. 
a higher best-fit temperature by 0.03 keV, which is larger 
than the statistical error quoted in Tabled However, the 
ratio of O to Ne was still found to be solar or larger (see 
Sect. 5). Inclusion of ROSAT All-Sky Survey (RASS) data 
into the analysis of the background allowed us to constrain 
better the Local Bubble component, e.g. a temperature of 
0.10 ± 0.01 keV was found but with a somewhat lower lu- 
minosity. Using these values instead of the standard ones 
quoted above had no effect on our results. 

In Fig.^we show the location of the spectral extracting 
areas. We select regions ^ 40' from the Coma center to the 
North- West, North, North-East, South and South-East. 
The South- West direction is complicated by the presence 
of the infalling subcluster (NGC4839), while some other 
pointings were affected by background flares, which after 
screening lead to insufficient exposures. 

3. Results 

Our major results are shown in Fig.|21and listed in Table 
We find that the soft excess in Coma on spatial scales of 
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30' — 50' is characterized by thermal plasma emission, be- 
cause of the clear detection of O lines, with a characteristic 
temperature of ^ 0.2 keV. Bonamente et al. (2003) found 
the same temperature characterizing the ROSAT PSPC 
observation of the Coma soft excess. The highest tem- 
perature of this component corresponds to the hot spot in 
Coma found in ASCA observations (Donnelly et al. 1999), 
and identified as an accreting zone of the Coma cluster. 



We have performed several internal checks of our re- 
sults. First, Fig. 121 shows the confirmation of our findings 
for the Coma- 11 field by analysis of the XMM-Newton 
MOSl data. The medium filter was also used for the MOS 
during these observations. For this analysis we have done a 
simple background analysis that uses the fiux in the MOSl 
field for the same instrument settings and background con- 
ditions but a different Nh- Since the MOS temperature of 
the hot component is consistent with that determined by 
the pn, we fix the temperature of that component to its 
pn value in the MOS spectral analysis. The MOS fit then 
agrees with the pn results for the intensities of the hot and 
warm components and for the temperature and element 
abundance for the warm component. 



Second, to entirely exclude the effect of the instrumen- 
tal or particle background, we have attempted an in-field 
background subtraction for the pn. Using the difference 
in the spatial distribution of the Coma hot and detec- 
tor background components within the Coma- 11 field, we 
were able to remove both of them using data only from the 
Coma-11 field. In particular, the region south of the point- 
ing center contains more of the cluster emission relative to 
the detector background and vice- versa for a region north 
of the pointing center. Since the spectrum of the CXB is 
similar to that of Coma hot emission (see Fig.EJ, it could 
be removed by over-subtracting the Coma hot flux. We 
exclude energies corresponding to the detector emission 
lines from the comparison, as their distribution is differ- 
ent from that of the continuum (e.g. Lumb et al. 2002). 
The resulting spectrum is displayed in Fig.[21and could be 
fitted with a single MEKAL model with 0.22 ± 0.01 keV 
temperature and a factor of 3 higher intensity than the 
Milky Way halo plus Local Bubble emission. Third, we 
have tried several compilations of the XMM background 
(Lumb et al. 2002; Read 2003) and found that the soft 
component in the Coma emission is clearly seen using any 
of them. 
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Fig. 2. Detailed decomposition of the pn spectrum (small grey 
crosses) in the Coma- 11 field into foreground and background 
components, obtained from the analysis of the APM field, plus 
hot and warm emission from the Coma cluster. Large crosses 
are the result of the in-field subtraction of Coma emission and 
detector background, possible due to differences in the spatial 
distributions between the Coma warm. Coma hot and detector 
background components. MW halo means Milky Way halo and 
CXB means cosmic X-ray background. 




0.6 0.8 1 2 

Energy, keV 

Fig. 3. MOS 1 APM field background subtracted spectrum of 
the Coma-11 field. Detailed modeling, accounting for the dif- 
ferences in Nh was not performed, but the difference is small, 
given the lower sensitivity of MOS CCDs at low energies com- 
pared to the pn. 

At the moment, it is difficult from the O line redshift 
alone to decide whether this emission is Galactic or ex- 
tragalactic. However, some apparent differences with the 
Galactic emission are already noticeable. The soft com- 
ponent is centered on Coma and has an amplitude ex- 
ceeding the variation of the underlying Galactic emission 
(Bonamente et al. 2003). An O abundance of 0.1 times 



solar is much lower than the Galactic value of one solar 
(Markevitch et al. 2002; Preyberg & Breitschwerdt 2002). 
The temperature variations of the Coma warm component 
are not seen in the quiet zones of Galactic emission at high 
latitudes. The observed warm component exceeds the level 
of the Milky Way halo plus Local Bubble continuum (at 
0.57 keV) emission by a factor of 3. Also the strength of 
the observed O lines are a factor of three larger than the 
corresponding values for the Galactic emission, and to sig- 
nificantly change the derived O abundance of the Coma 
warm emission a possible variation in the metallicity for 
the Galactic emission should be by a factor of ~ 3, which 
is not observed (Freyberg & Breitschwerdt 2002). As we 
will demonstrate below, ROSAT PSPC data strongly con- 
strain the contribution of the low-temperature Galactic 
components to the observed O line flux. 

~i ' 1 ' 1 — ' — I — ' — I 
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Fig. 4. ROSAT All-Sky Survey spectrum of Coma-11. Grey 
crosses show the background spectrum 2° north of the Coma 

center. Black crosses show the background subtracted spec- 
trum extracted at the position of the Coma-11 field. Lines 
show the best fit spectrum from the analysis of pn data folded 
through the ROSAT PSPC with no adjustment to the pn best 
fit. Dashed line is the warm component, dotted line is the hot 
component and the solid line is the sum. The best fit model to 
the pn data provides a statistically acceptable fit to the RASS 
data (x^ = 9.8 per 9 d.o.f.). 

4. Comparison with the ROSAT All-Sky Survey 

We use the RASS to provide an external check of our 
results. We note, however, that the spatial and spectral 
resolution and energy band of the ROSAT detectors are 
significantly different from those on XMM so care needs 
to be taken when performing this check. 

We first compare the fluxes implied by our analy- 
sis to those in the RASS data as presented in Snowden 
et al. (1997). We choose the 0.5-0.7 keV band (ROSAT 
R4 band), corresponding to position of the O lines and 
the 1.3-2.0 keV band (R7), where the relative contribu- 
tion of the hot component is the largest. We make the 
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comparison for two fields: Coma-11 and a field 4° North 
of Coma, which we call Sky Background. Several effects 
should be taken into account to provide a proper compar- 
ison between the XMM-Newton and ROSAT data. First, 
the X-ray mirrors of XMM-Newton have a 7% efficiency 
to stray-light (Lumb et al. 2002) with uncertainty in the 
energy dependence of 2%. This results in higher fluxes 
seen by the XMM detectors near bright objects and may 
lead to object-specific fiux variations. We use the RASS 
data to estimate the effect of stray-light pollution in the 
background and Coma fields by assuming uniform az- 
imuthal and radial dependence of the stray-light coming 
from XMM off-axis angles of 0.4° - 1.4°. This compo- 
nent should be subtracted from the XMM fluxes in the 
comparison. Second, the RASS maps of the diffuse X-ray 
background have point sources removed, which reduces 
the intensity of the background. We use the limits of 
Snowden et al. (1997) and the logN — logS of Hasinger 
et al. (2001) to estimate this removed flux. This compo- 
nent should be added to the RASS flux in the comparison. 
Third, RASS data on the sky background on the scales of 
the XMM field of view have rather large statistical uncer- 
tainties and some averaging should be carried out in order 
to provide a meaningful comparison. Fourth, the R4 im- 
ages exhibit variations larger than the statistical errors, 
so we quote the errors corresponding to the systematic 
variations. Fifth, in the analysis of XMM spectra a sepa- 
rate power-law component was introduced to remove the 
increased detector background due to soft protons (De 
Luca & Molendi 2002). Studies of this component by sev- 
eral groups (Katayama et al. 2002; Snowden S., Molendi 
S. 2003, private communications) show that its intensity 
could be determined by the spectrum at high energies (10- 
15 keV) where the mirrors have no collecting power for X- 
rays, while the instrumental background originates from 
energy losses of particles on the detector. From the com- 
parison of the XMM predicted and the RASS observed 
intensities of the Sky Background presented in Table |21 it 
is clear that our background model is applicable for the 
Coma cluster. Our results for Coma-11 are also in agree- 
ment with the RASS flux, although the later has a rather 
large uncertainty. 

The pixel size in the Snowden maps is rather large 
(12'), so we have performed a second analysis by extract- 
ing a spectrum from the raw RASS data at the exact lo- 
cation of the Coma-11 field, a circle of 8' radius centered 
on a 194.66521°, 5 = 28.406924°. We have extracted a 
background spectrum using a 0.6° x 2° region (long axis 
perpendicular to radius vector from the cluster center) 2° 
North of the cluster. We display the results in Fig. ^ We 
find that our pn model provides an acceptable fit to the 
RASS data with no adjustable parameters (x^ = 9.8 per 
9 degrees of freedom) . 

Fig. 21 also shows that a simple comparison of the coun- 
trates between the two instruments is not straightforward, 
particularly at energies below 0.5 keV. The ROSAT back- 
ground is very high in this energy range compared to the 
pn. This is caused by two differences between the instru- 



Table 2. Comparison between XMM-Newton Predictions of 
the ROSAT Surface Brightness from Observations (Coma-11) 
or Background Model (Sky Back) and Observed ROSAT All- 
Sky Survey Data. 



Field XMM-Newton pn RASS 
name predicted stray-light data removed flux 
lO""* PSPC-C counts s"^ arcmin"^ 



0.5-0.7 keV 






Coma-11 2.7 ±0.1 0.10 


2.0 ±0.5 


0.2 


Sky Back 0.66 ± 0.08 0.05 


0.5 ±0.1 


0.2 


1.3-2.0 keV 






Coma-11 3.2 ±0.1 0.14 


3.2 ±0.7 


0.2 


Sky Back 0.56 ± 0.04 0.05 


0.4 ±0.02 


0.2 



ments. At an energy of 0.25 keV the effective area of the 
ROSAT PSPC is about a factor of two higher than the 
XMM pn (~ 175 cm^^ vs. ^ 80 cm^^) while the energy 
resolution of the PSPC is about a factor of two worse 
than the pn (~ 1 vs. ^ 2) (Snowden et al. 1997 for 
the PSPC, XMM Users Handbook for the pn). At 0.25 
keV all components are dominated by the Local Bubble. 
The poor energy resolution of the ROSAT PSPC redis- 
tributes the larger number of counts from this component 
to higher energies. On the other hand, since the XMM- 
Newton Coma observations were made with the medium 
filter, most of the counts at lower energies are actually re- 
distributed from higher energies. This is why our results 
are not so sensitive to the model of the Local Bubble. The 
combination of these two effects produces substantial dif- 
ferences in the ratio of the Coma counts to the background 
counts below 0.5 keV in the two instruments and a proper 
comparison requires detailed spectral modeling as Fig. 0] 
shows. 

The comparison of Fig.|5]and Fig. ^reinforces our pre- 
vious analysis using the Snowden maps, but with a finer 
spatial and spectral resolution. For example, the inten- 
sity of the total background continuum at 0.57 keV in 
both figures is about equal to the Coma warm continuum. 
However, the unlimited field of view of the RASS allows 
a direct measurement of the background in Fig. 01 com- 
pared to a model of the background in Fig. [3 Thus we 
conclude that an incorrect background model producing 
the observed Coma Warm component is ruled out. 

Our analysis of the ROSAT data can also be used to 
set limits on the neutral gas fraction in a filament. From 
Table 1, the HII column is about 10^^ cm~^. From Fig. El 
the limit on the column of HI above the Galactic value 
is < 10^° cm~^. Hence, the HI fraction is less than 10%. 
Marginal indication of such absorption is present in Fig. 01 
as the data points are slightly below the model prediction. 

5. Interpretation 

In the rest of this Paper we present an interpretation as- 
suming an extragalactic origin of the warm component. 
With the CCD-type spectrometry we can place limits on 
the possible redshifts of the warm component, which for 
the Coma-11 field is determined as 0.007 ± 0.004 ± 0.015 
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(best-fit, statistical and systematic error). Any extra- 
galactic interpretation should therefore concentrate on the 
large-scale structure in front of the Coma cluster since the 
allowed redshift is 0.000 - 0.026 compared to the Coma 
value of 0.023. We have carried out an analysis of the CfA2 
galaxy catalog (Huchra et al. 1995) towards the Coma 
cluster, excising the South- West quadrant, where strong 
influence of the infalling subcluster NGC4839 on the veloc- 
ity distribution has been suggested (CoUes & Dunn 1996). 
Fig. illustrates our result, indicating a significant galaxy 
concentration in front of the Coma cluster, with velocities 
lying in the 4500 — 6000 km s~^ range. The sky density of 
the infall zone is shown in the insert of Fig.|5| It was con- 
structed from the number of galaxies in the 4500 — 6000 
km s~^ velocity range minus the number of galaxies in the 
8500 — 10000 km s~^ bin, corresponding to a similar dif- 
ference in the velocity dispersion from the cluster mean. 
By excising the galaxies at velocities lower than 6000 km 
s~^, we recover a projected velocity dispersion of ~ 750 
km s~^. Using the M — a relation of Finoguenov et al. 
(2001), we find that such velocity dispersion is more in 
accordance with mass of the Coma cluster. The excess of 
galaxies exhibits a flat behavior within 0.8 degrees (1.3 
Mpc), followed by a decline by a factor of 5 within 1.7 
degrees (2.8 Mpc) and a subsequent drop by two orders 
of magnitude within 10 degrees (16.7 Mpc) The soft emis- 
sion from the Coma is detected to a 2.6 Mpc distance 
from the center (Bonamente et al. 2003) in a remarkable 
correspondence to the galaxy filament. 

In the absence of the finger-of-god effect (decoupling 
of galaxies from the Hubble flow) , validated below by our 
conclusion on the filamentary origin of this galaxy con- 
centration, the infall zone is characterized by a 20 Mpc 
projected length. This length assumption affects the den- 
sity estimates, while the O abundance is only based on 
the assumption of coUisional equilibrium. We have verified 
the later assumption by studying the ionization curves for 
OVII and OVIII presented in Mathur et al. (2003). We 
have concluded that the assumption of pure coUisional 
ionization is valid for our data, since rig > 10~^ cm^'^ 
and T > 2 X 10^ K (0.17 keV). An advantage of our 
measurement is that we also determine the temperature 
by the continuum. Lower temperatures, which at densi- 
ties near 10^^ cm^"^ result in similar line ratios for OVII 
and OVIII, fail to produce the observed (H-e and He-e) 
bremsstrahlung flux at 0.7-1 keV. 

An important question we want to answer is whether 
the detected emission originates from a group or from a fil- 
ament. We cannot decide from the electron density of the 
structure (~ 50(/Xe//ip)/baryonPcrit), as it suits both (for 
the low metallicity of filament we take /Xg = 1.1, fj.p = 1.2 
and /baryon = 0.16 to corrcspoud to the first WMAP re- 
sults in Spergel et al. 2003). However, the implied mass of 
the structure is comparable to the mass of the Coma clus- 
ter (see also Bonamente et al. 2003). On the other hand 
the temperature of the emission is ^ 0.2 keV, almost two 
orders of magnitude lower than that of the Coma cluster. 
It takes a few hundred groups with virial temperature of 



0.2 keV to make up a mass of the structure, which leads 
to an overlapping virial radii if we are to fit them into the 
given volume of the structure. Thus we conclude that the 
observed structure is a filament. 



T 1 1 1 1 1 1 1 r 




Fig. 5. Galaxy distribution in the direction to the Coma clus- 
ter from CfA2 survey. Central 2 degrees in radius field centered 
on NGC4889 are shown, excluding the 60 degrees cone in the 
South- West direction from the center, to avoid the effect of the 
NGC4839 subcluster. An excess of galaxies over the Gaussian 
approximation of the velocity dispersion of the Coma cluster is 
seen in the 4500 — 6000 km s~^ velocity range. The insert shows 
the spatial distribution of the excess galaxies in the 4500 — 6000 
km s~^ velocity bin over the galaxies in the 8500 — 10000 km 
s~^ velocity bin, this time out to 10 degrees. 

For the Coma-11 field, where the statistics are the 
highest, we investigated the effect of relaxing the assump- 
tion of solar abundance ratios. We note that an assump- 
tion for C abundance is important for overall fitting and a 
solar C value would affect our results. When left free, how- 
ever, the C abundance tends to go to zero. Also, there is 
no systematic dependence of our results on the assumed 
C abundance as long as the C/0 ratio is solar or less, 
which is the correct assumption from the point of view of 
chemical enrichment schemes and observations of metal- 
poor stars in our Galaxy. Significant abundance measure- 
ments are: 0/0© = 0.14 ± 0.02, Ne/Ne© = 0.14 ± 0.06, 
Fe/Fe© = 0.04tg °i^ (assuming Fbq/H = 4.26 x lO'^ by 
number), indicating that Fe is underabundant by a fac- 
tor of three in respect to the solar Fe/0 ratio, imply- 
ing a dominant contribution of SN II to Fe enrichment. 
Element abundances for the hot emission obtained from 
our XMM data on the center of Coma are Mg/Mg© = 
0.36 ± 0.12, Si/Si© = 0.45 ± 0.08, S/S© = 0.01 ± 0.15, 
Fe/Fe© = 0.20 ± 0.03, also suggesting prevalence of SN 
II in Fe enrichment, although Fe abundance is a factor 
of 5 higher compared to the filament. The filament Ne/0 
ratio reveals a subtle difference with the OVI absorbers. 
Nicastro et al. (2002) reports Ne/0= 2 times solar for 
their X-ray absorption observations of material associated 
with the OVI absorbers. We observe a lower (solar) ratio 
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at 95% confidence. Lower ratios seem to be a characteris- 
tic of enrichment at ^ > 2 (Finoguenov et al. 2003b), while 
simulations suggest that the WHIM originates at z < 1 
(Cen & Ostriker 1999). 

While we would rather have more observational evi- 
dence on the dispersion of O to alpha element ratios at 
different sites, we believe that this is a major source of 
systematics in interpreting the element abundance of X- 
ray filaments as a universal value. To illustrate the point, 
we calculate the density of baryons traced by the OVI 
absorbers under two sets of assumptions. Scaling the orig- 
inal value of nhiOVI) = 0.0043/ifo^ of Tripp ct al. (2000) 
for ionization equilibrium implied by measurements of 
Mathur et al. (2003) and using the measurements of the 
O abundance reported here yields: 



Of- = 0.020/.- (^) 



Ifjovi)]- 

32 



(1) 



The formal errorbar is 0.006, mostly from the uncertainty 
in the estimate for ionization. If, on the other hand, our 
measurements of Ne abundance are used with the Ne/0 
ratio for OVI absorbers from Nicastro et al. (2002) then, 



0.14 lOl^'^/^h /([/(Oy/)]-i) 



IQlNe/H] 



32 



(2 



Therefore the second set of assumptions must be invalid 
since the total baryon density is fijj°'^' = 0.039, leaving no 
room for other major components of local baryons such 
as Lja absorbers (0.012 ± 0.002) and stars and clusters 
of galaxies (~ 0.006; e.g., Finoguenov ct al. 2003b and 
references therein). O depiction onto dust grains in the 
OVI absorbers, as suggested by Nicastro et al. (2002) as 
an explanation of the high Nc/0 ratio, docs not explain 
the unacceptably high baryon abundance implied by Eq. 
(2), since the solar Ne/0 ratio in our observations may be 
simply explained by dust sputtering. 

The observational determination of scaling relations 
between X-ray properties, such as luminosity Lx- gas tem- 
perature T, and entropy is crucial in establishing the phys- 
ical properties of the ICM. The slopes of the Lx — T 
and mass-temperature relations (e.g. Markevitch 1998; 
Finoguenov et al. 2001) and the gas entropy level (e.g. 
Ponman et al. 1999; Finoguenov et al. 2002) are all at vari- 
ance with model predictions based on pure gravitational 
heating (Kaiser 1986) and require the introduction of ex- 
tra physics to describe the thermodynamics of the ICM 
(e.g. Evrard & Henry 1991; Kaiser 1991). An often dis- 
cussed piece of extra physics is preheating of the baryons 
before they accrete onto the cluster. Most baryons that 
accrete onto clusters are thought to come from filaments, 
so we now have an opportunity to compare the thermody- 
namic properties of the filament gas with that of group and 
cluster gas. The entropy of the X-ray emitting filament 
is 150 keV cm^, which could be reproduced by heating 
while falling onto a filament (e.g. Cen & Ostriker 1999). 
In fact, estimating the expected Mach number of the ac- 
cretion shock when the filament enters the Coma cluster 



yields a value of 4, similar to predictions of simulations by 
Miniati et al. (2000). Regardless of the origin of the en- 
tropy of the filamentary gas, its entropy is much smaller 
than the 400 keV cm^ implied by ASCA observations of 
the outskirts of groups (Finoguenov et al. 2002), ruling 
out a universal preheating value. We note that the energy 
of the gas (|A:T = 0.36 ± 0.02 keV/particle) is similar, 
although higher than the SNe energy associated with en- 
richment of gas to the observed O abimdance (0.22 ± 0.03 
keV/particle), so heating by galactic winds is not ruled 
out. 

Although, the temperature of the filament is presently 
a factor of 40 lower than the Coma virial temperature, 
as it falls into the cluster it will be shock heated, and its 
final adiabat will be higher than in the self-similar case 
(Dos Santos & Dore 2002; Ponman et al. 2003). Since the 
mass of the filament is comparable to that of the Coma 
cluster, the combined object will also deviate from cluster 
scaling relations. Observations indicate that the gas- to- 
dark matter distribution will not be affected (Sanderson 
et al. 2003), but the temperature is higher for a given 
mass (Finoguenov et al. 2001). This process is probably 
not universal, since some groups of galaxies have shallow 
sgas profiles, thus possibly pointing to adiabatic compres- 
sion rather than shock heating of the preheated gas, as 
proposed by Tozzi & Norman (2001). 

If the structures reported here are only associated with 
large clusters of galaxies, their contribution to the baryon 
budget is negligible (0.1-1%). Of greater importance is 
to investigate the accreting environment of the massive 
groups and poor clusters that have a significant entry 
in the baryon budget (6%). However, early enrichment 
epoch, suggested by low Ne/0 ratio, makes X-ray fila- 
ments a substantial entry in the metal budget at high 
redshifts. 

The thermodynamic state of the filamentary gas causes 
a global feedback effect on the embedded galaxies by 
strangling the gas accretion (Finoguenov et al. 2003a: Oh 
& Benson 2003). The resulting star- formation will pro- 
ceed by consumption of the previously accumulated gas, 
in either quiescent mode as in the disk, or merger-induced 
bursts leading to formation of the spheroid (Somerville 
et al. 2001). Galaxy mergers are frequent inside galaxy 
groups (Kodama et al. 2001), but in a filament the infall 
of the gas will primarily be recorded in the star-formation 
history of the disk (Kennicutt et al. 1994). A relevant ob- 
servation could therefore shed light on the feedback epoch, 
which is crucial in understanding the relation between 
the X-ray filaments and OVI absorbers. The recent Sloan 
Digital Sky Survey (SDSS) discovery of passive spirals, in 
the same filament in front of Coma (Goto et al. 2003), is 
exactly what is expected from this strangulation process. 
The existence of these passive spirals lends further sup- 
port to the association of the soft X-ray excess with the 
Coma filament. As passive spirals are starting to be found 
in the outskirts of many clusters, the presence of the soft 
X-ray excess might be a wide-spread phenomenon among 
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the massive clusters, as indicated by first results on XMM- 
Newton REFLEX-DXL survey (Zhang et al. 2003). 

6. Conclusions 

We have discovered oxygen line emission associated with 
the outskirts of the Coma cluster. These data show that 
the previously observed soft X-ray excess in this cluster 
is due to thermal emission from a warm gas, possibly 
the long-sought WHIM. Typical parameters characteriz- 
ing the emitting material are a temperature ^ 0.2 keV, 
abundance ~ 0.1 solar, and density ~ 50/baryonPcriticai- 
The thermodynamical state of the gas in the Coma fil- 
ament reduces the star-formation rate of the embedded 
spiral galaxies, providing an explanation to the presence 
of the passive spirals in the SDSS survey in this region 
(Goto et al. 2003). 

Acknowledgements. The paper is based on observations ob- 
tained with XMM-Newton, an ESA science mission witli in- 
struments and contributions directly funded by ESA Member 
States and the USA (NASA). The XMM-Newton project 
is supported by the Bundesministerium fur Bildung und 
Forschung/Deutsches Zentrum fiir Luft- und Raumfahrt 
(BMFT/DLR), the Max-Planck Society and the Ifeidenhain- 
Stiftung, and also by PPARC, CEA, ONES, and ASI. AF ac- 
knowledges receiving the Max-Plank-Gesellschaft Fellowship. 
Authors thank Konrad Dennerl for applying the experimental 
bias correction software to the Coma-11 observation. AF prof- 
ited from discussions with Michael Freyberg, Margaret Geller, 
Fabrizio Nicastro, Maxim Markevitch, Monique Arnaud, Doris 
Neumann, Alexey Vikhlinin, and Francesco Miniati. AF would 
like to particularly mention a continuous support of this work 
by Richard Lieu and his collaborators, Max Bonamente, John 
Mittaz, and Jelle Kaastra. JPH thanks MPE for its hospitality 
during several extended visits. We acknowledge an anonymous 
referee for pointing out an importance of the detailed compar- 
ison between XMM-Newton and RASS. 

References 

Arnaud, M., Aghanim, N., Gastaud, R., et al. 2001, A&A, 
365, L67 

Bonamente, M., Lieu, R., Joy, M., & Nevalainen, J., 2002, 
ApJ, 576, 688 

Bonamente, M., Joy, M., & Lieu, R. 2003, ApJ, 585, 722 
Briel, U.G., Henry, J.P., & Bohringer, H. 1992, A&A, 259, 
L31 

Briel, U.G., & Henry, J.P. 1995, A&A, 302, 9 
Briel, U.G., Henry, J.P., Lumb, D.H., et al. 2001, A&A, 
365, L60 

Cen, R., & Ostriker, J.P. 1999, ApJ, 514, 1 
Colles, M., & Dunn, A.M. 1996, ApJ, 458, 435 
Dave, R., Cen, R., Ostriker, J.P., et al. 2001, ApJ, 552, 
473 

De Luca, A., & Molendi, S. 2002, preprint 

[astro-ph/0202480 
Donnelly, R.H., Markevitch, M., Forman, W., ct al. 1999, 

ApJ, 513, 690 
Dos Santos, S., & Dore, O. 2002, A&A, 383, 450 



Evrard, A.E., & Henry, J.P. 1991, ApJ, 383, 95 

Fang, T., MarshaU, H.L., Lee, J.C., Davis, D.S., & 

Canizares, C.R. 2002, ApJ, 572, L127 
Fang, T., Sembach, K. R., & Canizares, C.R. 2003, ApJ, 

586, L49 

Finoguenov, A., Reiprich, T., & Bohringer, H. 2001, A&A, 
368, 749 

Finoguenov, A., Jones, C, Bohringer, H., & Ponman, T.J. 

2002, ApJ, 578, 74 
Finoguenov, A., Borgani, S., Tornatore, L., & Bohringer, 

H. 2003a, A&A, 398, L35 
Finoguenov, A., Burkert, A., & Bohringer, H. 2003b, ApJ, 

594, in press, (astro-ph/0305190) 
Freyberg, M. J., & Breitschwerdt, D. 2002, Proc. of 

JENAM 2002, in press 
Fukugita, M., Hogan, C.J., & Peebles, P.J.E. 1998, ApJ, 

503, 518 

Goto, T., Okamura, S., Bernardi, M, et al. 2003, PASJ 
subm., (astro-ph/0301303) 

Haberl, F., Briel, U.G., Dennerl, K., & Zavlin V.E. 2002, 
to appear in Proceedings of 'New Visions of the X- 
ray Universe in the XMM-Newton and Chandra Era', 
ESTEC, The Netherlands; astro-ph/0203235 

Hasinger, G., Altieri, B., Arnaud, M., et al. 2001, A&A, 
365, 45 

Huchra, J.P., Geller, M.J., & Corwin, H.G.Jr. 1995, ApJS, 
99,391 

Jansen, F., Lumb, D., Altieri, B., et al. 2001, A&A, 365, 
LI 

Kaastra, J. S., Lieu, R., Tamura, T., Paerels, F.B.S., & 

den Herder, J. W. 2003, A&A, 397, 445 
Kaiser, N. 1986, MNRAS, 222, 323 
Kaiser, N. 1991, ApJ, 383, 104 

Katayama, H., Takahashi, I., Ikebe, Y., Matsushita, 
K., Freyberg, M.J. 2002, A&A submitted 
(astro-ph/0210135) 

Kennicutt, R., Tamblyn, P., & Congdon, C. 1994, ApJ, 
435, 22 

Kodama, T., Small, L, Nakata, F., Okamura, S., Bower, 

R.G. 2001, ApJ, 562, L9 
KuU, A., & Bohringer, H. 1999, A&A, 341, 23 
Lieu, R., Mittaz, J.P.D., Bowyer, S., Breen, J.O., 

Lockman, F.J., Murphy, E.M., Hwang, C.-Y. 1996, 

Science, 274, 1335. 
Lumb, D.H., Warwick, R.S., Page, M., & De Luca, A. 

2002, A&A, 389, 93 
Lumb, D.H., Finoguenov, A., Saxton, R., ct al. 2003, 

SPIE, 4851, 255 
Markevitch, M. 1998, ApJ, 504, 27 

Markevitch, M., Bautz, M.W., Biller, B., et al. 2003, ApJ, 
583, 70 

Mathur, S., Weinberg, D.H., & Chen, X. 2003, ApJ, 582, 
82 

Miniati, F., Ryu, D., Kang, H., et al. 2000, ApJ, 542, 608 
Nagamine, K., Fukugita, M., Cen, R., & Ostriker, J.P. 

2001, ApJ, 558, 497 
Neumann, D. M., Arnaud, M., Gastaud, R., et al. 2001, 

A&A, 365, L74 



Finoguenov et al.: XMM-Newton discovery of an X-ray filament in Coma. 



9 



Nicastro, F., Zezas, A., Drake, J., et al. 2002, ApJ, 573, 
157 

Oh, S.P., & Benson, A.J. 2003, MNRAS 342 664 
Ponman, T.J., Cannon, D.B., & Navarro, J.F. 1999, 

Nature, 397, 135 
Ponman, T.J., Sanderson, A.J.R., & Finoguenov, A. 2003, 

MNRAS, 343, 331 
Rasmussen, A., Kahn, S. M., & Paerels, F. 2003, ASSL 

Conf. Proc, eds. Rosemberg J.L., & Putnam, M.E., 

Kluwer, 281, 109 
Read, A.M. 2003, preprint |astro-ph/ 02124 36 
Sanderson, A.J.R., Ponman, T.J., Finoguenov, A., Lloyd- 

Davies, E. J., & Markevitch, M. 2003, MNRAS, 340, 

989 

Scharf, C, Donahue, M., Voit, G.M., Rosati, P., & 

Postman, M. 2000, ApJ, 528, L73 
Smith, R.K., Brickhouse, N.S., Liedahl, D.A.,& Raymond, 

J.C. 2001, ApJ, 556, L91 
Snowden, S.L., Egger, R., Freyberg, M.J., et al. 1997, ApJ, 

485, 125 

Somerville, R.S., Primack, J.R., & Faber, S.M. 2001, 

MNRAS, 320, 504 
Striider, L., Briel, U.G., Dennerl, K., et al. 2001, A&A, 

365, L18 

Tozzi, P., & Norman, C. 2001, ApJ, 546, 63 
Tripp, T.M., Savage, B.D., & Jenkins, E.B. 2000, ApJ, 
534, LI 

Wang, Q. D., Connolly, A. & Brunner, R. 1997, ApJ, 487, 
L13y 

Zappacosta, L., Mannucci, F., Maiolino, R., et al. 2002, 
A&A, 394, 7 

Zhang, Y.-Y., Fi noguenov, A., Bohr inger, H., et al. 2003, 
A&A subm., (|astro-ph/0306115|| 



